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SUMMARY 


A method for computing the turbulent heat transfer to two- 
dimensional and a xi symmetric bodies in high-speed flow is derived from 
the integrated equations of the boundary layer by utilizing a Stewart son 
type transformation. A modified Reynolds analogy between skin friction 
and heat transfer which depends upon local pressure gradient results 
from the analysis. Exact and approximate solutions are derived from the 
approximate differential equations; the exact solution is applicable for 
arbitrary initial (transition) conditions and the approximate solution 
requires fully developed turbulent flow from stagnation point or leading 
edge . 


For the exact solution, the assumption of fully developed turbulent 
flow from the stagnation point yields a solution representative of a 
variety of arbitrary transition solutions. The exact solution (restricted 
to stagnation initial conditions) and the approximate solution are shown 
to agree within 5 percent when applied to several blunt shapes. The 
present solutions generally predict the measured heating rates on these 
bodies within the accuracy of the measurements except where transition 
was thought to begin in the region of the peak predicted heat transfer. 

The present solutions appear to be sufficiently accurate for design 
purposes. The exact solution alone offers the generality of arbitrary 
transition; the approximate solution offers accuracy and relative sim- 
plicity but requires the assumption of stagnation -point initial condi- 
tions. Because of the scarcity of knowledge of the transition phenomenon, 
this latter restriction does not seriously impair the usefulness of the 
approximation . 
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INTRODUCTION 


Reentry conditions at high Reynolds number and the possibility of 
surfaces roughened by collisions with meteoric particles and dust are 
conducive to the existence of turbulent boundary layers on missiles, 
satellites, spaceships, and similar vehicles entering a planetary 
atmosphere. A primary concern of the ballistic missile program, for 
example, is the accurate prediction of turbulent heat transfer to a 
blunt body in high-speed flow. This problem has inspired a great deal 
of research in recent years; although many approximate theories have 
evolved, none appears to be sufficiently general for widespread 
application. 

The problem is the computation of turtulent heat transfer from a 
compressible (possibly dissociated) gas to a cold wall in the presence 
of a pressure gradient. Except for the case with dissociation, the 
zero-pressure-gradient skin friction may be computed with sufficient 
accuracy either from mixing length theories, for example, those of 
Van Driest for a plate (ref. l) and for a cone (ref. 2), or from an 
empirical power law such as the Blasius law. (See, for example, ref. 3.) 
Compressibility in the latter case may be taken into account through a 
reference temperature formulation of the power law. In both methods, 
heat transfer is computed from the skin friction through the use of 
Reynolds analogy between heat transfer and skin friction suitably modi- 
fied for the effect of a nonunit Prandtl number. 

The presence of a pressure gradient seriously modifies the flow and 
the subsequent friction and heat-transfer behavior. Most recent methods 
of predicting the turbulent heat transfer to a body in the presence of a 
pressure gradient center about extensions cf the inverse power law friction 
coefficient. Some of these are described by Libby and Cresci in refer- 
ence 4. The simplest uses flat-plate skin friction and heat-transfer rela- 
tions evaluated for the local flow conditions on the body in question. 

Rose, Probstein, and Adams (ref. 5) combine a flat-plate shear law with 
the complete momentum equation integrated across the boundary layer and 
including pressure gradient. Heat transfer is computed from the resulting 
skin friction through the flat -plate Reynolds analogy modified to account 
for Prandtl number. Reference 5 reasons that the effect of the pressure 
gradient on the velocity and temperature profiles is small and thus the 
flat-plate friction law and analogy between heat transfer and skin fric- 
tion are valid. 

A relation between heat transfer and stin friction including the 
effect of pressure gradient may be derived from the integrated momentum 
and energy equations. (See refs. 4, 6, and 7.) Bloom and Martellucci 
(ref. 6) were apparently the first to use this relation, and they combined 
it with an inverse power friction law and a reference enthalpy correction 
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for compressibility effects. The results of this method and a similar 
method (ref. 4) indicated predicted heat transfer lower than that found 
experimentally, and excessive dependence of the modified Reynolds anal- 
ogy upon pressure gradient. Assumption of the validity of the Crocco 
form of the total enthalpy -velocity relation in the boundary layer 
appears to have caused this trend. 


Another method of accounting for compressibility and pressure gra- 
dient was developed by Reshotko and Tucker. (See ref. 8.) This method 
utilizes the Stewartson transformation (ref. 9) applied to the turbulent 
case, a friction law including the effect of pressure gradient upon 
shape factor, and a Reynolds analogy factor evaluated from an approximate 
solution for laminar flow in reference 10. This approach is restricted 
to isothermal walls with adverse and small favorable pressure gradients 
regardless of wall temperature and to flow with large favorable pressure 
gradients over highly cooled surfaces. M. Richard Dennison while at 
Missile Systems Division of Lockheed Aircraft Corporation (in a paper 
not generally available) extended Van Driest' s analysis to include the 
effect of pressure gradient upon Reynolds analogy, neglecting, however, 
the direct effect of pressure gradient upon the momentum equation. This 
assumption appears to be justified only for highly cooled surfaces. The 
analysis is further restricted to isothermal walls. 

The present analysis attempts to correct the deficiencies of the 
aforementioned methods and provides a comprehensive method for the com- 
putation of turbulent heat transfer in compressible flow for high sub- 
sonic and supersonic as well as hypersonic conditions. The effect of 
pressure gradient is retained in the integrated momentum equation and a 
modified Reynolds analogy is derived. The Stewartson transformation is 
used to reduce these equations to a form close to the incompressible 
relations. An incompressible friction relation which neglects the 
effect of pressure gradient upon the boundary-layer shape factor is 
assumed to be valid in the transformed coordinate system and the com- 
pressible friction coefficient is then calculated by using a form of 
the reference temperature method. The heat transfer is then calculated 
from the friction and the modified Reynolds analogy. 

This treatment includes the possibility of nonisothermal heated or 
cooled surfaces but is limited to favorable and small adverse pressure 
gradients . Although dissociation is neglected, variable specific heats 
may be approximated through use of a linear interpolation formula. 

Heat-transfer rates predicted by the present method are compared 
with experimental results and with other predictions. 
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Subscripts: 
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enthalpy convection thickness, / —(l - £v 
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transformed, enthalpy convection thickness, / -H- / l _ 
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evaluated at wall in insulated case 
evaluated outside "boundary layer 

evaluated at reference location (stagnation point) outside 
boundary layer 

evaluated at reference location in boundary layer 
evaluated at transition point 
evaluated at wall 
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Primed quantities are evaluated at the temperature T 1 . A bar over 
a quantity represents evaluation in the transformed plane; for example, 
C f* V ^ bar under a quantity means it has been made dimension- 
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THEORY 


Derivation of the Transformed Equations 

The equations of motion and energy in the physical plane .- The 
boundary -layer equations of motion and energy may be integrated across 
the boundary layer to obtain the classical Karman momentum integral and 
its energy counterpart. See, for example, references 5 and 7 . The 
resulting equations are for axisymmetric steady flow; 
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Equations (l) and (2) are very general within the usual Prandtl boundary- 
layer assumptions requiring that the boundary -layer thickness be small 
compared with the body radius r and that the fluid be composed of a 
single constituent gas of arbitrary properties. Air may be assumed to 
fulfill the last-named requirement if the concentration of its constitu- 
ents are fixed (no dissociation) . 


The momentum and energy equations have been solved for many cases 
of laminar flows. Insufficient knowledge of the form of the turbulent 
shear , heat transfer, and boundary- layer profiles has made analytic 
solutions cf these equations impossible for turbulent flow. With the 
aid of certain empirical relations , however, approximate solutions to 
the momentum equation have been obtained. (See, for example, refs. 1, 

3, and 3.) From these solutions it is possible to obtain local friction 
and heat-transfer coefficients, the latter being computed through use 


of the Reynolds analogy Ng^ 


2 


which is valid for zero -pres sure - 


gradient flow with unit Prandtl number and constant wall temperature. 
For a nonunit Prandtl number, the relation if altered to the form 

Ng t = f^Np r ) This problem has been treated most recently in 


reference 11. 
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A generalized relation between heat transfer and skin friction, 
hereafter called the modified Reynolds analogy, may be obtained by com- 
bination of equations (l) and ( 2 ) with the result 

_ q * - 2 21 + d ( He * 9) _ cpfi + 

P e u e H e * 0 2 H e * dx \ 0 /u e dx 

A similar modified Reynolds analogy was first derived in reference 6 and 
has been used in references k and 7- Equation (k) reduces to the con- 
ventional Reynolds analogy for zero pressure gradient with an isothermal 
wall and equal to 1. 



Application of the Stewartson transformation .- Mager has shown 
(ref. 12) that the Stewartson transformation can be applied to turbulent 
flow with unit Prandtl number and zero heat transfer in order to reduce 
the equation of motion to its incompressible form. Where the Prandtl 
number differs from 1 and heat transfer is permitted, the transformation 
reduces the momentum equation to a form close to but differing from that 
for incompressible flow. In the present analysis, a modified form of 
the Stewartson transformation is applied to the momentum equation and 
the modified Reynolds analogy relation. 


The transformation used herein is given by the following relations: 


X(x) 



p r p r 
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Vy Jo 


dy 
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The subscript r refers to conditions at some local reference 
within the boundary layer. The subscript o refers to conditions at 
some arbitrary reference location outside the boundary layer. Also, 

P r H r = f(x) 

p Q p 0 = Constant 
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From the definition of the thickness functions (9 and qp) and the 
relation between the transformed and physical normal coordinate (eq. (6)) 
the following relations may be obtained: 

p o/ h o\ 1//2 - 

6 = — — 0 (8) 

PeWI K 
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In a like manner, the displacement thickness 6* is written as 


6 * > 


SI Rc 


o 1 - f f 


-/* 

'-'O 


1 - 1 dY 


(11a) 


In order to account approximately for the effects of variable spe- 
cific heats, the following linear relation between density and enthalpy 
is assumed to be valid: 
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where a is a constant adjusted so that equation (l 2 ) is exact at the 
wall. Fay and Riddell (ref. 13 ) used a correlation formula similar to 

equation (12) but including, in addition, a term in ^1 - ^ • The 

coefficient of this added term was small compared with a, and that term 
is herein neglected. 


Combination of equations ( 8 ), (ll), (12), and the definition of 
total enthalpy 

[u(x,y)] 2 


H(x,y) = h(x,y) + 


(13) 


yields 
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(14a) 



The momentum equation (eq. (l)) and the modified Reynolds analogy 
(eq. (4)) may he transformed to the Stewartson plan through the use of 
equations ( 5 ) to (l4). If the energy equation outside the boundary 
layer is 
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the momentum equation is 
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(16) 
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and the modified Reynolds analogy is 
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Equation (l6) takes on a form almost like the incompressible form 


with the definition of the transformed friction coefficient c*. as 


c„ p (I p 
f o o e 

2 p r p r p Q 


( 18 ) 


The transformed friction coefficient is related to the transformed wall 
shear t w by the definition 


= T w 

2 “ p Q u e 2 


(19) 


Combination of equations (l8), (19), (6), and the definition of the 
friction coefficient in the physical plane = T v/ P e u e 2 ) ylel^ 8 the 

following relation between wall shears 
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( 20 ) 


In an analogous manner, the transformed heat -transfer coefficient is 
defined as 


% _ q w p e D o^o 

Po^e^e* Pe u e^e* Po ^rMr 


The heat transfers in the two systems are reLated by 


^w 


PoPo/ h o\ ly/2 
PrPr\h e j ^ 


( 21 ) 


( 22 ) 


Substitution of definitions (l8) and (21) in:o the transformed momentum 
equation and modified Reynolds analogy (eqs. (l6) and ( 17 )), respectively) 
reduces these relations to their final general form. 
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The momentum equation is 
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and the modified Reynolds analogy is 

% 


© c f - 

= J- + CD 


?o U eV 9 2 


+ <P < 


1 ffl e + 1 d(<p/0) 


H g * dX 9/0 dX 


1 + — + 


2 Hg 


1 + — (a - l) 


a 


V 62 

H P b 


J_ fUe 

U e dX 


(21+) 


Equation (23), even for constant specific heat (a = l) , is not in 
general equivalent to the incompressible momentum equation because of 


the presence of the term -a. 


5 He U e dX ) 


For unit Prandtl number 


with zero heat transfer, H e * vanishes; thus, the coefficient of the 
velocity gradient term would take on its incompressible form 2 + — 
in that case. 


The present analysis for axisymmetric flows can be specialized to 
the two-dimensional case by considering r as a constant and dr/dx = 0; 
thus, r is eliminated from the momentum equation. 

Relation between total enthalpy and velocity profiles . - In order 
to solve the system of equations outlined in the preceding section, it 
is necessary to specify velocity and total enthalpy profiles and a fric- 
tion law, but it is desirable to investigate first a relation between 
total enthalpy and velocity profiles which should have rather general 
significance . 

For the special case of zero pressure gradient, unit Prandtl number, 
and an isothermal wall, the Crocco relation between total enthalpy and 
velocity is valid. This equation is 


H = a + bu 
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Boundary conditions at the wall and stream give the specific relation 

H* _ _u_ 

H e * ' u e 

Pressure gradient, nonunit Prandtl number, and nonisothermal wall will 
each probably alter this relation. In reference 11 , for example, with 
the assumption of nonunit Prandtl number in the laminar sublayer, the 
total enthalpy is found to be a quadratic function of velocity in the 
sublayer. 

Dennison assumed the validity of the quadratic relation 

H = a + bu + cu^ (25) 

for the case with nonzero pressure gradient but with unit Prandtl number 
and isothermal wall. The coefficients of the velocity terms were found 
to be functions of the local pressure gradient and the heat transfer 
was evaluated from the skin friction by usir.g a relation analogous to 
the modified Reynolds analogy of the present report. A multiplication 
factor was then suggested to account for nor unit Prandtl number. 

Rather than complicate the present analysis by combining the Prandtl 
number and pressure gradient effects into tie enthalpy-velocity relation, 
the approach of the present report is to follow Dennison in assuming 
that the Prandtl number is 1 and that equation (25) is valid through the 
entire boundary layer. The coefficients of the velocity terms in equa- 
tion (25) are then evaluated as functions of pressure gradient and wall 
enthalpy. A correction to account for nonurit Prandtl number is later 
introduced. Boundary conditions for equaticn (25) were such that at 
the wall 

y = Y = 0 


u = U = 0 


H = IV 


and at the outer edge of the boundary layer 

y = 6 

Y = A 

u = u e 

U = U e 
H = H e 
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The form of equation (25) and the second boundary condition require 
equal velocity and enthalpy thicknesses, that is, 5 U = 6 ^ = 6 and 

= ^ = A. This restriction is not serious since a unit Prandtl 
number has already been assumed; thus equal or nearly equal thicknesses 
are implied. Application of the boundary conditions to equation (25) 
yields 



(26a) 


= — fl 




(26b) 


Equations (26) are assumed to be valid independent of the type of 
velocity profiles and the friction lav describing the flov. A qualita- 
tive check of their validity may be obtained by comparison with laminar- 
flow results. Plotted in figure 1 are results of the similar solutions 
of reference 14 for zero and favorable pressure gradients and cooled 
walls. The quantity (3 of reference 14 is a pressure gradient parameter 
which is zero for zero pressure gradient and increasingly positive for 
increasingly favorable pressure gradients. Plotted also is equation (26a, 

XI 2 

for various values of c -S— from 0 to 1.0, values which also represent 

H e* 

zero and favorable pressure gradients and cooled walls. The quantities 
u 2 

6 and c — — are not simply related to one another and the present com- 

V 

parison is intended as only a qualitative one. With this qualification 
the two sets of curves in figure 1 show sufficient resemblance in shape 
to lend support to the present approach. 

With the aid of equation (26b) the thickness functions and cp 

may be expressed in terms of integrals of velocity functions as follows: 
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Definition of Z as 


He* 9 


yields the relations 


h = h + z i 

00 ^ 


(30) 


l + z 


(31) 


One more equation is needed to complete the formulation of the 
problem. Equation (24) represents the relation between heat transfer 
and skin friction in the transformed plane. Another relation between 
the same two quantities, one which must be compatible with equation (24), 
may be derived from equation ( 26 b) as follovs: 


Consider the derivative of enthalpy with respect to velocity from 
equation (26b) evaluated at the wall. This relation is 





(32) 


The shear and heat transfer at the wall can, in general, be written 
respectively, as 
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W 



(33a) 
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Y=0 


(331) 


because the apparent turbulent shear and heat transfer vanish at the 
wall. The values of p., k, and c are not here specified but apply 

ir 

to a common location. The ratio of shear to heat transfer is easily 
obtained from equations (33) as 
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and, since the Prandtl number has been assumed to be equal to 1, 
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Combining equations (32) and (34) yields 


P 0 U e H e* 


i . zip: 
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(35) 


Equations (24) and (35 ) > each a form of the modified Reynolds anal- 
og y, are assumed to be equivalent. A differential equation for Z 
results from this assumption. After substitution of equations ( 30 ) 
and ( 31 ) into the result, the following equation is obtained 
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A s r - 1 + a 5s- Z-2. = 1 + 


_e_Z±=l + ^^ + (a - 1 ) 
H e X H e 0 
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Note that r is a function of Z but that A is not. Finally, by 
using equation ( 30 ) the momentum equation ( 23 ) is rewritten as 
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d0 + / r_ + 1 
dX lu e dX r 

Equations (36) and (38) form a pair of simultaneous nonlinear differ- 
ential equations of the first order. The boundary conditions chosen are 
those resulting from the simple assumption of instantaneous transition 
from laminar to fully developed turbulent flow; that is, 0 = 0^, 

Z = Z t at X = X t . 

The special case of constant specific heats yields the simpler 
equations 



_ 5sl zl 

H e 9 Hg \ 


( 39 a) 


K. 8 -, 

A = 1 + — — (39b) 

Hg 0 

1 dH e* 

The terms in have been retained in equation (36) in 

Hg* dX 

order that the restriction to isothermal walls may be relaxed to include 
H\# =>■»<*>• The assumption is made that the gross effects of a non- 
isothermal wall are accounted for in this eqaation without the necessity 
of further modifying the enthalpy-velocity relation (26a) or the depend- 
ent relation between q w and Cf (eq. (35)). 


Application of a Simple Velocity Profile 
and Friction Law 

Assumption of power velocity profile ani friction law .- In the 
absence of detailed knowledge of the turbule it velocity profiles and 
friction law for compressible flow about a b idy of revolution with 
streamwise pressure gradient, certain assumptions are made. These are: 

(1) The effect of compressibility is completely accounted for by 
the Stewartson transformation; thus, incompressible flow relations may 
be assumed to be valid in the transformed plane. 

( 2 ) The effect of pressure gradient upon the velocity profile and 
friction law is negligible; thus, zero-pressure -gradient relations may 
be used. The incompressible, zero -pres sure -gradient relations chosen 
for the present analysis are 
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u m 1 /" 

(40) 

u e ' W 

C„ fo u ‘eX -1 / m 

4 = 

(41) 


where A^ is assumed to be constant. Velocity profiles of the form of 
equation (40) have been observed in flows with pressure gradient. (See 
ref. 15 .) The Ludwieg-Tillman friction law of this reference is of the 
form of equation (4l) but with A m as a function of shape factor. 

Since the shape factor is not altered seriously by favorable pressure 
gradients but is a strong function of adverse pressure gradients, the 
present analysis is limited to favorable and weak adverse pressure 
gradients . 


( 3 ) The velocity profile and friction law (eqs. (40) and (4l), 
respectively) are unaffected by the curvature in the axially symmetric 
case so that these two-dimensional relations may be assumed to be valid. 
With the assumption that the boundary-layer thickness is small compared 
with the body radius, the local curvature due to the axial symmetry 
would be expected to have a small influence upon the local friction law 
because, locally, the turbulent exchange should far outweigh the momentum 
change caused by the varying cross-sectional area. On these physical 
grounds the assumption of no effect of body radius upon profiles and 
shear law seems to be justified. Some methematical justification might 
be seen from the following: If a turbulent analog to the laminar Mangier 

transformation were employed, equations (40) and (4l) would be assumed 
to be valid in the resulting two-dimensional flow. It can be shown 
that such a procedure leads to an expression for the physical (three- 
dimensional) friction coefficient equivalent to that obtained by the 
present method, that is, the expression which results from combination 
of equations ( 5 ) , (6), (l8), and (4l), 




P U 

■cr o e 


(42) 


Thickness functions .- For the velocity profile of equation (40), 
the thickness functions of equations (lOa), ( 11a) , and (29) may be 
evaluated 
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(45) 


Differential equations .- Equation (4l) is substituted into the 
momentum equation (38) and, after some rearranging, the result is 
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a(e) m , 

dX 


/ i ■ nH-1 

v m / lu e dX r dX 



(44) 


The modified Reynolds analogy equation (3 6) becomes 
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(45) 


The dependent variables are now (0) 


tt\ m 


anc. 


Z. 


Nondimensionalization . - Equations (44) and (45) may be made dimen- 
sionless with respect to the constant reference conditions p Q , p Q; 

p , h , and a constant reference length Z. Dimensionless quantities 
are then, for example, 
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dx 


Equations (44) and (45) become 
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(47) 


where 



(48a) 

(48b) 


The quantity © is a "reduced" momentum thickness, and the Reynolds 
number dependence of the equations is implicit in its definition and the 
Reynolds number enters these equations in no other way. Equations (46) 
and (47) are nonlinear differential equations through the dependence 
of the quantity r upon the dependent variable Z. Step-by-step numer- 
ical solution is required and the solution is hereafter termed the 
"exact " solution . 
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Evaluation of Skin Friction ana Heat Transfer 


When the momentum thickness parameter © and the factor Z have 
been computed, the skin friction and heat transfer may be obtained as 
follows: Combination of equations (4l) and (48) gives 



(49) 


In the physical plane, from equations (l 8 ) and (4-9), 
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(50) 


Then, for unit Prandtl number, from equatior.s (21) and ( 35 ) 
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(51) 


For Prandtl number other than one, the enthalpy driving potential for 
heat transfer is H aw - H w , and this substitution for Hg* in equa- 
tion ( 51 ), along with the empirical Prandtl number correction, gives 
the result 


where 


1 

Nstfao )"* 1 



P M* 

-r-r 
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-e 
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ITH-1 




-q, 
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St p u H * 
H e e aw 


(52) 


Special Cases 

Solution for H e * equal zero.- FOr the case H e * equal to zero 
with arbitrary pressure gradient and constar.t specific heats, equa- 
tions (46) and (47) may be solved in closed form. The term — — 

H e* d£ 
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in equation (77) appears to be indeterminate but may be seen to be zero 
since for any constant value this term is identically zero 

regardless of how small the difference H e - H w = H e * is allowed to 
become . 


The differential equations describing this flow are, from equa- 
tions (37), (46), and (47), 


d© ( m + 1 \ f r d ^e 1 _ (m + 1 

d| \ m )\U d^ r d 5/ m ( m 
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(54) 


r = 2 + 


a = 1 + -± 

0 


(55a) 

(55“b) 


On q 

The thickness ratios — and ~ were evaluated for power profiles 

0 A 

and are given in equations (73) • Since n was assumed to be constant, 
these ratios, P, and A are then constant, and equations (53) and (5*0 
are linear and may be integrated in closed form. The momentum equa- 
tion (53) may be written after some manipulation as 
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Equation (56) is integrated to give 
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24 


In the special case 0^ = = 0, 


© = A 


m + l\/T, \ m JO 


I mfi r(itH-l)-! 

(r) m (Uj m d| 


tth -1 r(ntf l) "1 
(X) m W " 


The modified Reynolds analogy equation ( 54 ) is of the form 


d£ 1 Ld? 


and the solution is 


Z = e 




where 


Mi) = 




i + ±\ - ^2 

I Ho d l 


In general, the integral f f 1 (|)d| will have to be evaluated numeri- 
cally after completion of the solution of 1 he momentum equation. For 
the special case © t = | t = 0; however, the- solution for © yields 


ro+1 r(m+l^ -1 

m (g) m (^e) m = __m_ d jgga g 

m + 1 ft m+1 r(m+: ) -1 m + 1 d£ 

(r) m ( 14 ) 111 d| 
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and 


n . ffi ( 1 + J_ N 

/ f l(i )d ^ = (p l) nH ' 1 \ A 




-(r-i) 


( 59 ) 


Equation (59) is substituted into equation (58) and the following equa- 
tion is obtained: 


Z = 
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d| + K 


(60) 


The constant K must be set equal to zero in order that equation (60) 
leads to a stagnation-region solution consistent with that developed in 
the appendix. The limits on the integral are set from zero to g, which, 
in conjunction with K = 0, allows Z to approach zero as expected far 

dlL, 

downstream on a blunted cylinder or cone as -2H. approaches zero. 

d| 


The final result is 


Z 



m 

m+1 


(F) 




i 


0 



(61) 


Equations (5Tb) and (6l) provide explicit solutions to the momentum and 
modified Reynolds analogy equations for the case 0^ = = 0 for 

H e * = 0, which, for unit Prandtl number and recovery factor, corresponds 
to an insulated wall. 

Other special cases .- The special cases of the flat plate, the 
cone in supersonic flow, and two- and three-dimensional stagnation 
regions are discussed in the appendix. The stagnation-region solutions 
have value in establishing initial conditions for the complete solution 
of equations (46) and (47) when the flow is assumed to be turbulent from 
the stagnation point. The results for two- and three-dimensional stag- 
nation regions are given in the appendix (eqs. (A10) to (A15) ) and plots 
of the variation of Z Q with the cooling ratio are shown in 

figures (2a) and (2b) for a = 1, m = 4, and n = 7. Although the 
results (eqs. (AlO) to (AI5) ) are also valid for heated walls H^/Hg > 1, 
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only the cooled-wall results are shown. For stagnation flow, the reduced 

m-1 

momentum thickness © varies directly with | m , or in the physical 

plane, 0 varies directly with x 1 ’* . The momentum thickness at the 
stagnation point is zero, a consequence of the assumed relation between 
c f and 0. The term Z is independent of x in the stagnation region. 


Approximate Solution 

The special case of flow with pressure gradient but zero enthalpy 
difference H e * was shown to give simple closed-form solutions for 

e t = It = 0 and constant specific heats. ( See eqs. (57b) and (6l).) 

An approximate method for the general case bised on this special case 
is now developed. 

The nonlinearity of equations ( 46 ) and (47) is caused by the depend- 
ence of T upon Z. (See eq. ( 37 a )-) Figure 3 shows plots of the 
variation of T with cooling ratio H^/^e for various values of Z 

with a = 1 and n = 7. The upper limit of Z of 0.8 was chosen since, 
for this value, the Stanton number is zero for a finite friction coeffi- 
cient. (See eq. (52).) A value of Z in excess of this maximum would 
imply heat transfer in a direction opposite ;o that of the heat-transfer 
potential a behavior not expected on physical grounds except per- 

haps where H aw * passes through zero. Figure 3 shows that r is not 
a strong function of Z for given Hy^Hg. For Hy|H e approaching 
zero, the dependence is a maximum; for unit YL^j H e , the dependence dis- 
appears. Similar behavior would occur for hot walls (h w /H 6 > 1.0 
nonunit a. is not expected to alter this dependence seriously. 

The approximation assumes that T is given by equation ( 37 a ) spe- 
cialized for an average value of Z equal to 0 . 4 . Thus, 



H** , x u e 2 / 

z -°- k s = 2 + % j * (a - 11 sc 1 + i 


!4!!i 

He 0 


Equation (62) is assumed to be correct for ar arbitrary value of n. 
The momentum equation (46) is linear with F substituted for r and 

jji dU 

may be integrated upon the neglect of with respect to The 
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derivative 


dr 

d| 


is seen to depend upon 


dK/He) 

d| 


and 


da 

Hp 


both of which 


are assumed to be small. Because approaches zero downstream on 

blunted cylinders or cones, the approximation may fail in these regions. 
However, with this assumption and some manipulation, equation (46) is 
rewritten as 
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(63) 


which is identical to equation ( 56 ) except for the substitution of T 
for P. Solutions are therefore given by equations (57 a ) or (57^) f° r 
appropriate initial conditions with r replaced by 7. 

The energy equation (47) is solved in the same manner. When T 
is replaced by r, equation (47) becomes a linear differential equation 
of the form 


where now 


f x (6) 
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( 65 b) 


The general solution is 


Z = e 


P 


ft i(s)ae 


fo(5)d| + K 


(66) 


As in the special case, equation ( 66 ) must be evaluated numerically for 
arbitrary transition because of the form of the function © (eq. ( 57 a ))« 
For the special case of 0^ = = 0, however, a simple analytic solu- 

tion may again be derived. Thus, when the method used in the case 
H e * =0 is followed here, the result is 
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In the interest of simplicity, the case of arbitrary transition 
(eqs. (57a) and (66)) for the approximate solution will not be consid- 
ered further. The approximate solution is given by equations (57b) 
and (67) and is therefore restricted to the initial conditions 

- 0 , = Z Q or fully developed trrbulent flow from the 

leading edge. 

Stagnation-region solutions for the approximate case may be obtained 
from equations (57b) and (67). For a three-dimensional stagnation 
region, 
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and for a two-dimensional stagnation region. 
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The three- and two-dimensional stagnation values of Z Q as given by 

equations (69) and (71), respectively, have also been plotted on fig- 
ures 2(a) and 2(b) , for a = 1, n = 7; and m = 4 and do not differ 
significantly from the appropriate (negative root) values for the 
"exact" solution evaluated in the appendix. 

For clarification of the two methods discussed herein, a summary 
is presented in table I. 


Reference Conditions 

The constant reference conditions p Q , \± Q have been used as the 

pertinent properties for the flow in the transformed (incompressible) 
plane. This "incompressible" reference is assigned to the stagnation 
point in the case of a blunt body and local stagnation conditions for 
bodies with sharp leading edge. This requirement specifies isentropic 
flow along the body streamline in the inviscid flow for the latter case. 

The reference conditions in the boundary layer p r , p r are some- 
what more arbitrary. The original transformation (ref. 9) for laminar 
flow used the relation pp = 1, whereas reference Ik used the wall con- 
ditions p^,^. The choice of reference conditions in the laminar case 

influences the form of the transformed equations but not the solutions 
in the physical plane (except through the viscosity law chosen), since 
the solutions do not depend upon any assumptions as to the nature of 
the flow in the transformed plane. Because the turbulent solutions 
depend upon the assumption that incompressible relations are valid in 
the transformed plane, the reference conditions need to be carefully 
specified as follows. 

T w method.- If the laminar shear at the wall was considered and 
then transformed, this would give 

Tw * ^“(syjv 

= PwlV fb u\ 

p o p o h o °Ww 

If the quantity n ( — ] were interpreted as the transformed wall 

°W w 

shear r w , then equation (72) would be 
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( 73 ) 


Comparison of equations (73) and ( 20 ) indicates that the reference con- 
dition, for this interpretation, must be the* wall condition. For the 
flat plate this comparison gives, from equal. ion (A6), 

m n-1 

if = (74) 

c f \ p e^e/ \ P o/ 

T' method.- In reference l6, an intermediate temperature method 
for turbulent "flow is described. This method has been adapted to the 
flat -plate inverse power friction law as follows: Assume 


' w 


\ P'Ue 2 


= K n 


/p'u e x\ m " 1 


where is the coefficient of the equivalent incompressible law and 

p ! and p 1 are evaluated at an intermediate temperature T r . Then 
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The intermediate temperature is given for a .r by reference l6 as 


T 1 p . ' 

1- = 1 + 0.055M/ + 0.45 7 / - 1 
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\ e 


( 76 ) 
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The wall reference method (eq. (74)) is compared with the inter- 
mediate temperature method (eqs. ( 75 ) and ( 76 )) in figure (V) by means 
of plots of the variation of c fj^f with Mach number for wall-to- 

stream-temperature ratios of 0.2, 1, and 5 . For this figure it was 
assumed that the surface was at a temperature of 1,500° R for all Mach 
numbers, a fairly realistic condition for a reentry vehicle. Values of 
7 of 1.4 and m of 4 and the Sutherland viscosity law were used. For 
this case it is apparent that the two methods are significantly 
different . 


The same two methods are compared for the insulated wall case in 
figure 5 with a typical wind-tunnel condition imposed, namely, T Q 

equal to 540° R. A recovery factor equal to ( N p r ) was used with 

Np r equal to 0.72. Also shown on the figure are the experimental data 

of Coles (ref. 17 ) and Chapman and Kester (ref. 18) taken from refer- 
ence 16 . There is some tendency for the data to support the T f method 
over the T v method. Further support for use of the T 1 method is 

found in data recently obtained by Tendeland (ref. 19 ) for flow with 
heat transfer. 


The assumption is now made that the flat-plate T 1 method is cor- 
rect for arbitrary pressure gradient. The reference conditions for the 
present method are computed by requiring the equivalence of the T* 
compressibility correction to that of the present report (eq. (a 6 ) ) . 

The result of equating equation (a 6 ) with equation ( 75 ) is 


P r P r _ 

e 0 Po ” PqVW 


(77) 


with p T and \±' evaluated at T’ (eq. ( 76 )). However, for complete- 
ness, the T w method, given by equation ( 74 ), obtained from 


P r P r 

P M- 
o*o 


^w^w 

P P 
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(78) 


is evaluated as an alternate. 
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RESULTS AND DISCUSSION 


Method of Calculation 

Equations (46) and (47) for the "exact* method and equations ( 57b) 
and ( 67 ) for the approximate method were solved for various configura- 
tions on the IBM 704 card-programed computer at the Langley Aeronautical 
Laboratory. In order to simplify the programing, the equations were 
rewritten with the physical dimensionless coordinate x as the independ- 
ent variable as follows: 
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for the approximate case where 
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Because all available data were obtained in the moderate tempera- 
ture regime, the program was specialized to the case a = 1 (constant 
specific heats) and flow properties outside the boundary layer evaluated 
in the program from the pressure distribution by simple isentropic-f low 
relations with constant 7 . The viscosity was assumed to be given by 
the Sutherland law 

T 3/2 

M “ T + 198.70 R 

The program was arranged for arbitrary selected values of the constants 
A m , m, n, j&, 0/A, and Np^,. The enthalpy recovery factor r\ r 

was assumed to be equal to ^Np r ^ 1/3 

and the enthalpy difference H aw * 

was computed from 

H * = H e fh + T] u 2 _ It j 
aw eu-e > r -e / 


For purposes of comparison, both T* and T w methods for computing 
compressibility effects were included in the program. 

Unless otherwise specified in the subsequent discussion, the fol- 
lowing constants were used: 

Blasius law: 

m = k 

£ = 5 

x " 4 

7 = 1.4 N Pr = 0.72 


A m = 0.013 


l/ 7 th-power velocity profile: 

n = 7 


0 


9 

7 


Determination of Effect of Initial (Transition) Conditions 

In order to determine the effect of various transition conditions 
upon solutions of the "exact" equations, the discussion which follows 
assumed instantaneous transition from laminar flow at some location Xf 

For this discussion, the reference temperature as evaluated by the 
T’ method (eqs. (76) and (77)) is used. 
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No exact transition boundary conditions are known. One reasonable 

conditions is to assume that momentum is cor served: thus 9, 

7 laminar 

equals Q-turbulent at • An analogous procedure to determine Z^ 

is to conserve the enthalpy convection thickness cp. The purpose of 
figures 6 to 12 is to indicate the degree to which the results are 
affected by various transition boundary conditions. 


In figures 6, 7, and 8 are plotted Q/l, Z, and 




T - T ’ 
■ L w 1 aw 


respec- 


tively, as computed for a hemisphere in a F = 5 stream. The test is 
reported in reference 20 as run 57. Transition was assumed to occur at 




0.4, 0.6, and 0.8 with 



chosen as equal to the appropriate 


laminar value (computed from the similar solutions of ref. 14 by using 
a technique similar to that of Stine and Wanlass (ref. 21)). The term 
Z t was allowed to vary for each value of x t - For each x t the choice 

of Z t has little effect on the solutions for Z and ajll - T \ 

" / \ » aW ) 

except in the immediate vicinity of ^ . There is no noticeable effect 
upon 0 / Z at all. 


The choice of 


*t 


and 


(t: 


has greater influence upon the results ; 


ag may be seen in figures 9 to 11. Since the choice of Z t was shown 


to influence only weakly the solutions for given and 


a 


, for 


these three figures Z t was set equal to ths value of Z at z equals 
obtained from the solution for the stagnation initial conditions 

h = (V°’ 

given x , the choice of 


which are also shown in figures 6 to 12. For 


(i 


strongly affects the solutions locally, 
the effect diminishing downstream. In addition, increasing x+. increases 
the effect of varying that is, the effect is felt further down- 

stream, as expected. 

The data obtained in reference 20 appear to favor the solutions 
with large ^ with Zt fixed near x = 0.8 and to favor the solu- 
tions with (&\ = (& \ 


laminar 


at x^. =0.6 and 0.8 in the region 


X > 1.0. 
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The results of similar calculations for heat-transfer coefficient 
on a flat-faced body are shown in figure 12. The conditions are those 
of run 6l of reference 20. Transition was assumed to occur at = 0.2, 

0.4, and 0.6 for various values of . Because of the previously 

shown negligible effect of Z^., the values of Z^ were chosen in the 
same manner as for the hemisphere (figs. 9 to ll) • Figure 12 shows the 
same behavior for the flat-faced model as found for the hemisphere, 

that is, the effect of [^ ] decreases with increasing x and, when 

Wt 


= 0.2 or 0.4, does not appear to be important a moderate distance 
from the assumed transition point. The experimental data for this flat- 
faced body are plotted in figure 12 and appear to favor transition at 

equals 0.4 with equal to 0.00072 or 0 . 00108 . 


Some investigators have advocated starting the turbulent calcula- 
tion from the stagnation point with finite initial momentum thickness. 
(See, for example, ref. 5-) They reason that a nonzero boundary-layer 
thickness must exist at a stagnation point, such as the laminar value 
which is easily calculated. Solutions for the present method were 

obtained for a wide range of (j\ for x^ equal to zero for both 


spherical and flat -faced bodies. These results showed that the perturba- 
tions in Z, and - T &w ^ caused by the nonzero 

out extremely rapidly and are indistinguishable from those for - 0, 

x. = 0, downstream of x = 0.25- These results have not been shown in 


figures 9 to 12. Since the flow in the region close to the stagnation 
point is expected to be laminar, it must be concluded that for all 


practical purposes the assumption of finite 



at the stagnation 


point for the present method is practically equivalent to the assumption 


of zero 



Figures 6 to 12 show that the "exact" solution started at the stag- 
nation point fairly well represents the solutions begun at an assumed 
transition point, particularly for early transition. Therefore, this 
solution is useful when the location of transition is unknown. However, 
use of stagnation-point initial conditions does not imply that turbulent 
flow exists over the whole surface but merely assumes that the turbulent 
flow downstream of transition is approximately equivalent to a hypothet- 
ical turbulent flow which originates at the stagnation point. 
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Comparison of Present Method With Experiment 

Results of application of the present method are shown in figures 15 
to 19 corresponding to experimental data presented in references 20, 22, 
and 25 . Both exact and approximate solutiors were evaluated with the 
calculation restricted to the isothermal-wall case. This condition was 
very nearly met in all the experiments discussed. The values of H^/Hg 

chosen were approximate means of the local values and are indicated in 
table II, a summary of the experimental data used in the figures. In 
all cases the experimental pressure distributions presented in the refer- 
ences were used. 


The exact method was evaluated for stagnation initial conditions 
= = 0 > = Z o) ^ or one transition case in which x^. was 

chosen from the experimental data, ©*, was assumed equal to the laminar 
value and Z t was arbitrarily set equal to its value at x = x^ as 

given by the solution for stagnation initial conditions. Use of stag- 
nation initial conditions was implicit in the approximate solution. 

Plotted in these figures is the variation of the heat-transfer 

coefficient — - — with the dimensionless surface distance x = 

x w “ x aw l 

The reference length is the body radius for ^he hemisphere and flat- 
faced body of reference 20 (figs. 13 and 14) and the hemispheres of ref- 
erence 22 (figs. 15 and l6) , and the nose radius for the sphere cones of 
reference 23 (figs. 17 > 18, and 19 ) . In eacn case the (a) parts of the 
figures correspond to the T ? method, and the (b) parts of the figures 
to the T w method. ' 

Examination of the figures reveals the following: 

(1) In general, the results using the r ." method for c f jc f are 

slightly higher (generally less than 5 percent higher in the peak heating 
region) than those with c f j c f given by the T w method. The data show 

a slight preference for the T' method. 

(2) The approximate solution gives results within 5 percent of 
those of the exact solution using stagnation- point initial conditions. 
Neither of these methods predicts completely the experimental data, 
although most trends not too close to the transition region are satis- 
factorily predicted. 


( 5 ) The exact solution with transition initial conditions predicts 
the measured peak heat -transfer data within approximately 10 percent in 
those cases where transition appears to have begun upstream of the 
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theoretical peak (figs, 15, 18 , and 19) * In figures 15 and 18 , the 
location of the theoretical peak is slightly upstream of the experi- 
mental location. When transition appears to have begun close to the 
predicted peak heating region (figs. 13, l6, and 17)^ agreement at the 
peak may be either good (within 15 percent, fig. l6) or poor (figs. 13 
and 17) . In these three cases the peak heating data are thought to be 
transitional. 

( 4 ) Downstream of the peak heating regions agreement is generally 
within 10 percent except for the extreme downstream locations where 
heating rates are low and percent differences between theory and experi- 
ment tend to increase (figs. 13 and 15 to 19) . An exception is the flat- 
faced body where the predicted heating rates are 20 to 25 percent greater 
than the experimental rates over the entire face. 

{5) Using a transition initial condition in preference to stagna- 
tion initial conditions generally improves or leaves unchanged agreement 
between "exact" theory and experiment except in two of the cases where 
the peak heating data are thought to be transitional. (See figs. 13 
and 17. ) 

In regard to these comparisons, caution must be exercised in quan- 
titatively comparing data and theory as a 10 to 20 percent uncertainty 
is common in heat -transfer measurements. The question of choosing tran- 
sition conditions for the present method leads to further uncertainty. 
Although the technique of requiring conservation of momentum thickness 
may appear to be logical for an instability-induced transition, a 
roughness -induced transition might suggest a discontinuous momentum 
thickness. Figures 11 and 12 demonstrate how transition with 9 ^ 

larger than ^ laminar can ^- eac ^ to better agreement between theory and 

experiment for the same data as in figures 13 and 14 , respectively. In 
both of these cases, transition may be roughness induced. (See ref. 20 .) 

Comparison With Other Methods 

Comparison with Dennisons method . - The method of Dennison has 
been evaluated for the seven experimental configurations previously 
discussed. A simplified form which Dennison states is valid only on the 
forebodies of blunted cones and cylinders and invalid far downstream has 
been used along with a Prandtl number correction in the modified Reynolds 
analogy suggested by Dennison and used also in the present analysis. 

The results are shown in the (a) parts of figures 13 to 19 in each case, 
although Dennison’s method is not dependent upon a reference -temperature 
compressibility correction. The comparisons made in the following para- 
graph refer only to the present T r method and Dennison's method. 
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The method, of Dennison appears to predict heating rates about as 
well, in general, as does the present method with stagnation initial 
conditions. The present method with transition predicts the peak heating 
rates much more accurately than does Dennison’s method for the configu- 
rations of figures 15, l6, 18, and 19. For the hemisphere at M = 5 
(fig. 13) and the sphere cone (fig. 17), the Dennison predictions are 
closer to the peak heating rates than the present solutions, which over- 
predict the experimental peak heating rates However, these measurements 
appear to be in a transition region. Heating rates on the flat-faced 
body at M = 5 (fig- 14) are predicted by Dennison's method with excel- 
lent accuracy. 


Comparison with incompressible flat-plate formula . - In order to 
provide a familiar reference curve in figures 13 to 19, the heat-transfer 
rate was computed by using the Blasius incompressible flat-plate law and 
Reynolds analogy with the assumption of no effect of compressibility 
upon skin friction. The equations are 

tF - 0.0296(N Re )' I; 
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The results of using these simple relations are also plotted on the 
(a) parts of figures 13 to 19. This method however, is independent of 
compressibility. 

Although these results agree with the more complete results of the 
present analyses as well as with Dennison's method and data in some 
regions, differences are rather large over most of the surfaces. The 
occasional agreement is the result of comper. sating effects. The present 
method includes the effect of the favorable pressure gradient upon 
momentum thickness in the momentum equation. This tends to increase the 
skin-friction coefficient over the comparable flat-plate value. On the 
other hand, the effect of pressure gradient upon the Reynolds analogy 
factor is such that the ratio of heat-transfer rate to skin-friction 
coefficient is reduced below the comparable flat-plate value. In some 
regions these effects approximately cancel end agreement between the 
pressure gradient and flat -plate solutions Jesuits, particularly for the 
flat-nosed body. (See fig. 14.) In other legions one effect predomi- 
nates and the solutions do not agree at all. The compressibility cor- 
rection factor included in the present methcd further alters the compari- 
son. Because of this fortuitous part-time egreement, the incompressible 
flat-plate method is not suitable for blunt-nosed bodies. In the fig- 
ures discussed here, this method yields heat -transfer rates generally 
higher than the data or the results of the jressure-gradient methods. 
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The method of Bloom and Martellucci . - Bloom and Martellucci (ref. 6 ) 
use the integrated momentum and energy equations to derive a modified 
Reynolds analogy relation analogous to that of equation ( 4 ) of the pres- 
ent report. This relation is used in conjunction with an inverse power 
friction law, power velocity profile, reference -temperature compressibil- 
ity correction, and a relation between total enthalpy and velocity inde- 
pendent of pressure gradient in order to obtain the heat -transfer rate. 

The results of these calculations are not included in the present report 
but were shown in references 4 and 6 to underestimate seriously the 
turbulent heating rates on a blunt body. This behavior must be attrib- 
uted to the form of the enthalpy-velocity relation used. By not including 
the effects of pressure gradient in the total enthalpy profile, a modi- 
fied Reynolds analogy factor Pw/ T w which depends directly on velocity 
gradient is obtained, which is analogous to equation ( 4 ) herein. This 
relation shows that positive velocity gradients tend to reduce drasti- 
cally the Reynolds analogy factor and yield low heating rates. When the 
pressure gradient is included in the enthalpy profile, however, as in 
the present analysis, the modified Reynolds analogy factor depends 
instead upon an integral of the velocity gradient (for example, eq. (67) 
for Z) , and thus the effect of pressure gradient is reduced. This 
result is apparent in the generally satisfactory agreement between meas- 
ured heating rates and those predicted by the present method as shown 
in figures 13 to 19. 


Discussion 

Choice of initial conditions for the exact solution and validity 
of the approximate solution .- In a previous section it was shown that 
the choice of initial conditions at the transition point did not seri- 
ously affect the predicted heat-transfer rates a moderate distance down- 
stream of the assumed transition point. With no real knowledge of what 
transition conditions to use, figures 13 to 19 show that use of stagna- 
tion initial conditions can lead to reasonably accurate predictions of 
the turbulent heating rates a moderate distance downstream of the tran- 
sition point or region. It follows, then, that the approximate method 
presented herein may also be very useful, in spite of the fact that it 
requires stagnation initial conditions. Results of using this method 
also compare very favorably with the experimental data. 

If it is desired to utilize high-speed computing equipment in the 
calculation of turbulent heat transfer by the present method, the exact 
solution should be employed as it offers the generality of arbitrary 
initial (transition) conditions as well as exactness (within the frame- 
work of all the other assumptions) . Where no such equipment is avail- 
able, the approximate method offers a convenient alternative restricted 
to stagnation initial conditions. One might improve the accuracy of 
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the approximate method by a more judicious choice of a mean value of Z 
used in evaluating P (eq. (62)) but the need for doing so is not 
apparent in the present results. Caution must be used in applying this 
approximation, however, as it was derived by neglecting dZ/d£ with 
respect to dU^dg, , a questionable assumption in the essentially zero 

pres sure -gradient downstream regions of blunted cones or cylinders. In 
a subsequent section, this error is shown to be small for one case - a 
spherically blunted cylinder in shock-tube flow. 


An alternate friction law . - In order to determine the effect of 
changing the friction law, the heat-transfer coefficients for the hemi- 
sphere and flat -faced bodies at Mach number of 5 (ref. 20) were com- 
puted for the friction law 
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given in reference 3 and attributed to Falkner (ref. 24) . Results of 
these calculations are shown in figures 20 and 21 for hemispherical and 
flat-face bodies, respectively. Also shown are the results obtained by 
using A^j = 0.013 and m = 4 (modified Blasius law) as in the previous 

discussion, and the experimental results. Theoretical results are shown 
only for the special case where the flow is considered to be turbulent 
from the stagnation point = ©t = 0, Zt = Z 0 ^ for the exact solution 

(T‘ method). 


For a flat plate in incompressible flow, the modified Blasius law 
is considered to be a good approximation to the Prandtl-Schlichting 
logarithmic law for Reynolds numbers based on surface distances of from 
lO^ to 10?, and the Falkner law may be used to approximate the loga- 
rithmic law for Reynolds numbers of from 10^ to 10^-0. Approximate 
Reynolds number limitations such as these may be considered to be appro- 
priate for the calculation of heat transfer by the present method. For 
the experiments considered herein, Reynolds numbers based upon local 
external fluid properties and surface distance from the stagnation point 
never exceeded 2.5 x 10^; thus, the Blasius Law should be appropriate. 
Figure 20 for the hemisphere shows little difference in the two results, 
although the Blasius result is slightly closer to the data downstream 
of the transition region. Figure 21 for the flat-faced body shows that 
the Falkner result agrees with the data very closely, whereas the 
Blasius result averages about 20 percent high. This behavior remains 
unexplained. 

If it were desired to compute heating rates by the present method 
to a body for which the local Reynolds number based upon surface distance 
covered both of these ranges, the calculation could begin with the 
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Blasius law carried up to an appropriate location and then continue with 
the Falkner law for the remainder of the surface, with the requirement 
that 0 and Z were continuous at the "switchover point." The heat- 
transfer rate would he discontinuous at this point hut could easily he 
faired to obtain a turbulent heating rate continuous in x. 

Bloom and Martellucci (ref. 6) used an inverse power friction law 
with A m equal to O.OO636I and m equal to 5 * 92 , a law very similar 
to the Falkner relation. They also claimed to have computed one case by 
using a Blasius law (A m = 0 . 0128 , m = 4 ) and, when the results were 
compared, they found a difference of nearly a factor of 2 . They appear 
to have used a value of A m equal to 0.0225 in the actual calculations, 
however, and its use would account for the discrepancy. 

Validity of the approximate method in downstream regions . - The 
accuracy of the approximate solution in the downstream region of a 
spherically blunted cylinder was tested for a hypothetical shock-tube 
situation. The model was assumed to be one-half inch in diameter, and 
the traveling shock was assumed to move at a Mach number of 6.5 relative 
to the fixed wall and based upon the speed of sound of the undisturbed 
air, which was assumed at room temperature ( 540 ° F) and a pressure of 
40 centimeters mercury absolute. Local flow conditions around the model 
were found by assuming a modified Newtonian pressure distribution; 
stagnation conditions were computed from the charts of reference 25* 

For convenience in computing the heat transfer, the flow external to 
the boundary layer was assumed to be isentropic with 7 equal to 1.2. 
The wall temperature was assumed to be equal to room temperature, a 
valid assumption for the short -duration tests inherent in shock-tube 
use, and the computed stagnation temperature was 6,318° R. 

Figure 22 shows a plot of the variation of the ratio of local heat- 
transfer rate to stagnation heat -transfer rate (computed by the method 
of ref. 13) for exact and approximate solutions with x = x/Z, where Z 
is the body radius. Only the T r method is shown, and only stagnation 
initial conditions are chosen for the exact solution. Equation (76) 
for T 1 is assumed to be valid independent of 7. The region x < 1.57 
corresponds to the hemispherical nose, and the region x > 1.57 repre- 
sents the downstream cylinder. Figure 22 shows that the approximate 
solution differs from the exact by less than 5 percent over the entire 
surface considered, x ^ 3 . 0 . This result is surprising in the region 
where & > 1 . 57 , since in this region the Newtonian velocity gradient 
is zero whereas the gradient of Z is still finite. In obtaining the 
approximate solution, dZ/d£ was neglected as small compared with 
dU^dg, an assumption which is not true for x > 1.57- This relative 

insensitivity to the neglect of the dZ/d£ terms in the downstream 
regions somewhat widens the applicability of the approximate solution 
presented herein. 
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CONCLUDING REMARKS 


A method for computing the turbulent heat transfer to bodies in 
high-speed two-dimensional and axisymmetric flow fields was derived 
from the integrated equations of the boundary layer. The derivation 
utilized a Stewart son-type transformation of the equations, the assump- 
tion of a quadratic dependence of total enthalpy upon velocity across 
the boundary layer (in which the coefficients of the velocity terms are 
functions of pressure gradient), and a compressibility correction for 
skin-friction coefficient based upon flat-plate experimental data. The 
method was applied specifically for inverse -power-law velocity profiles 
and friction coefficients. Exact and approximate solutions were derived 
for these differential equations wherein the exact solution was valid 
for any initial conditions, and the approxinate solution required assump- 
tion of a fully developed turbulent flow frcm stagnation point or leading 
edge with appropriate initial conditions. A primary result of the 
derivation was the existence, within the present assumptions, of a modi- 
fied Reynolds analogy between heat transfer and skin friction for turbu- 
lent flow, a relation which depends upon strearnwise pressure gradient. 

Examination of the exact solution for various initial conditions 
(representing a variety of transition values) for flow about a hemisphere 
and a flat-faced body indicated that choice of the initial condition 
would influence the result close to the transition point but that this 
influence would decay in the downstream regions. It was shown that 
hypothetical initial conditions at a stagnation point or leading edge 
led to a solution which would satisfactorily approximate those solutions 
with arbitrary initial conditions downstream of their respective initial 
strearnwise location. 

Comparison of the present results with available experimental heat- 
transfer data on seven axisymmetric blunt bodies showed general agree- 
ment between experiment and theory, in most cases within the experimental 
accuracy. Use of a transition initial condition in preference to a 
stagnation initial condition tended to improve agreement between experi- 
ment and exact theory except for cases where transition appeared to 
begin in the region of peak predicted heat transfer. Results using the 
approximate solution never deviated from those of the exact solution 
with stagnation initial conditions by more than 5 percent. 

The present method appears to be capable of predicting turbulent 
heat transfer with sufficient accuracy for design purposes, except per- 
haps close to transition. Although the exact solution offers the gener- 
ality of arbitrary transition initial conditions, the computation is 
sufficiently involved to make desirable the employment of high-speed 
automatic computing equipment. In the absence of specific knowledge of 



the transition phenomenon, the assumption of fully developed turbulent 
flow from the stagnation point or leading edge still leads to highly 
useful predictions in the region where the flow is actually turbulent. 
With this assumption, the approximate solution has value in that it 
lends itself to hand calculation with little loss of accuracy. 

The present method is restricted to flow without dissociation. 
Variable specific heats may be included in an approximate manner; how- 
ever, no experimental data including this effect were available for 
comparison. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., September 30 , 1958- 



APPENDIX 


SPECIAL CASES 


Under certain conditions, equations (46' and (47) may be solved 
analytically. The case with H e * = 0 is treated in the text. Some 

additional examples are shown in the following sections. 


Flat Plate or Cylinder With Sharp Leading Edge, 

Constant Specific Heat 

For the flat plate or cylinder with shaj-p leading edge and constant 
specific heat, the equations are simplified by the restriction 
du^ dr 

— — - — =0, a = 1. Equation (46) is integrated immediately and the 
dx dx 

following equation results: 


0 = ®t + 





(Ala) 


If the flow is assumed to be turbulent from the leading edge, that is, 
®t = £ t = 0, then 


0 = 




(Alb) 


The skin friction for an all-turbulent flow is found from equations (4l) 
and (Alb) to be 


’i- 


(A2) 


One example of a friction law is the Biasius pipe-resistance formula as 
modified in reference 3 for which A^j = 0.013 and m = 4 and for which 
equation (A2) becomes 



1+5 



The constant K may be evaluated by consideration of the boundary con- 
ditions- When the simple case of constant wall enthalpy is considered, 


dH e* d (He - \) 

d| 


0 


and a particular solution of (Ah) is Z = 0. This is the solution of 
interest, the flat-plate Reynolds analogy. 

The relation between compressible and incompressible skin-friction 
coefficients is computed from the definition of X and Cf in conjunc- 
tion with equation (A2). It is assumed that reference conditions p r 
and y i will, in general, be functions of the local flow outside the 
boundary layer. For this flat -plate case, then, p r and p r are con- 
stant. If equations (5a) and (l8) are substituted into equation (A2) 
there results 
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The ratio of the skin frictions is, from equations (A 2) and (A5) , 

m 1 

Of _ ‘ “- 1 

Cf \ p e^o/ \^o ! 

at the same Reynolds number. 


(A5) 


(a6) 


Cone in Supersonic Flow With Attached Shock, 

Constant Specific Heats 

For the cone in supersonic flow with attached shock and constant 
specific heats, no pressure gradient exists for zero angle of attack; 

du e d ^ dr 

thus, = — — =0. In this case, however, rax and — = Constant. 

dx d£ dx 

The momentum equation (46) becomes 
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for which a solution, specialized immediately to 0^ = = 0, is 


© - 


0 =) 


m+1 
m d| 


(A7) 


Since rax, then r a £ because p r and p r are again constant < 
Equation (A7) becomes 


© = Jsl±jlVu x 


m 


(A8) 
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Combination of equations (4l) and (A8) yields 


£f 

2 



m + 



2m + 1 



1 

nH-1 


(A9) 


Just as in the case of the flat plate, the solution Z = 0 for 
constant wall enthalpy satisfies the differential equation for Z. 
Furthermore, the same compressibility correction for skin friction holds 
(eq. (A6)). Finally, it is apparent when comparing the two cases that 


( c l)cone _ ( C f) cone _ / 2m + A nH-l 

( Cf ) plate ( Cf ) Plate 


at the same Reynolds number and for the same wall and stream conditions. 
In other words, for the Same friction coefficient on cone and plate, 


w 


cone 


2m + 1 
m 


(%e) 


plate 


when wall and stream conditions are the sarnie on cone and plate. For the 
special case m = 4, this relation is 




plate 


a result which compares favorably with that of Van Driest (ref. 2) who 
obtained 


( N Re) cone 2 ( N Re) plate 


Stagnation Region on Two- and Three-Dimensional Bodies 

Although turbulent flow at a stagnation point has little physical 
significance, use of stagnation-point boundary conditions ^x^- = 0^ in 

the "exact" differential equations (46) and (47), as well as in the 
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approximate solution, necessitates a stagnation region solution. Such 
a solution is easily obtained for the exact case from equations (4b) 
and (47) by requiring that = a^, r = x for axisymmetric flow, and 

u = UqX for two-dimensional flow, where ag is the dimensionless 

st agnat i on -po int velocity gradient. By this procedure, for a three- 
dimensional stagnation region, 


0 = An 


a© 


( 6 A), 
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Z = Zq = 



(A10) 


(All) 


and for a two-dimensional stagnation region 


0 = Am 
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(A12) 




Two roots appear in each equation for Z Q . The negative root must 
be chosen as that which gives physically reasonable values of Z Q , that 

is, between zero and 0.8. (See approximate solution.) This is seen 
from figure 2, where Zq for the special case a = 1, n=7> m=U 

has been plotted against stagnation-point wall cooling ratio 

two- and three-dimensional stagnation points (eqs. (All) and (A15) ) • 
is seen that the negative root is appropriate. 

Equations (All) and (A13) also show that Z is not a function of | 
in the stagnation region, and hence = °* 
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TABLE II.- SUMMARY OF EXPERIMENTAL DATA USED TO COMPARE WITH 

PREDICTED HEATING RATES 


Figure 


13 


Reference 


Remarks 


Run 57; — = O.663; hemisphere; 
Hg 

Z = Radius 


Run 6l ; ^ = 0.660; flat -faced 
cylinder; Z = Radius of cylinder 

= 0.4-9; Rq = 19.26 x 10^ (nota- 
He 

tion of ref. 22); hemisphere; 

Z = Radius 


= 0.77; F d = 10.85 X 10 6 (nota- 
tion of ref. 22); hemisphere; 

Z = Radius 


Run 3-12; — = O.55; sphere-cone; 
Z = Nose radius 


Run 5-8; — = 0.60; sphere -cone; 
Z = Nose radius 


! 

23 

Hw . 

Run 7-3; — — = 0.60; sphere-cone; 
H e 

i 

Z = Nose radius 















Enthalpy difference ratio, H*/H 



Figure 1 .- Variation of enthalpy difference ratio with velocity ratio in 
boundary layer for present method and that of reference 14 . 





(a) Three -dimensional . 

Figure 2.- Variation of stagnation-point values of Reynolds analogy 
parameter Z q with wall cooling ratio = i; m = 4; 


Value of Reynolds analogy parameter at stagnation point, Zj (positive 
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Local Mach number, M e 

Figure 4.- Variation of fiat-plate skin-friction compressibility correc 
tion with local Mach number for various wall- to stream-temperature 
ratios. T w = 1,500° R. 
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Figure 5» — Variation of skin— friction compressibility correction with 
local Mach number for an insulated flat plate. T = 5b0° R. 
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Figure 10.- Distribution of Reynolds analogy parameter for hemisphere 
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Figure l6.- Comparison of theory and experiment for a hemisphere in a Mach number 2 stream. 











Experimental data (ref. 23) 

Exact solution, stagnation initial conditions 
Approximate solution 
Exact solution, Xt ~ 0* - 
Dennison 

Incompressible flat plate 
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(a) T' method. 


(b) T w method. 


Figure 19.- Comparison of theory and experiment for a sphere-cone in a 
shrouded tunnel. (Run 7-3 of ref. 23.) 
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Figure 20.- Comparison of results of using two friction laws with experimental data. Hemi- 
sphere in a Mach number 5 stream; T' method. 




Distribution of heat transfer 



Figure 22.- Distribution of heat transfer on hemisphere cylinder for 
typical shock-tube flow conditions. Shock Mach number equals 6.5 
pressure ahead of moving shock equals kO centimeters mercury; 

T ’ method . 
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